Charge and spin excitations of insulating lamellar copper oxides 
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A consistent description of low-energy charge and spin responses of the insulating Sr2Cu02Cl2 
lamellar system is found in the framework of a one-band Hubbard model which besides U includes 
hoppings up to 3'^'' nearest-neighbors. By combining mean-field calculations, exact diagonalization 
(ED) results, and Quantum Monte Carlo simulations (QMC), we analyze both charge and spin 
degrees of freedom responses as observed by optical conductivity, ARPES, Raman and inelastic 
neutron scattering experiments. Within this effective model, long-range hopping processes flatten 
the quasiparticle band around (0, tt). We calculate also the non-resonant Aig and Big Raman profiles 
and show that the latter is composed by two main features, which are attributed to 2- and 4-magnon 
scattering. 



I. INTRODUCTION 

Recent angular-resolved-photoemission spectroscopy 
(ARPES) measurements on the hard to dope insulat- 
ing Sr2Cu02Cl2 system have provided ,for the first time, 
data for the single-hole dispersion e(q) in an antiferro- 
magnetic background 0]. These data, as well as opti- 
cal absorption measurements , give information about 
charge excitations of the insulating cuprates. On the 
other hand, the spin excitations of the Cu02 planes have 
been tested by inelastic neutron and Raman scattering 
experiments [^,^. They show that the low energy spin 
excitations of insulating cuprates are well described by 
the two-dimensional spin-i antiferromagnetic Heisenberg 
(AFH) model Much theoretical work has considered 
both excitations separately, and therefore a simultaneous 
description of charge and spin degrees of freedom using 
the same model is lacking. 

Experimental results on the undoped Sr2Cu02Cl2 
lamellar cuprate 0-|4j, provide an unique opportunity 
to test, at the same time, the description of charge and 
spin responses as is obtained from current theoretical 
models for these strongly correlated systems In this 
work, we analize the electronic structure of insulating 
Sr2Cu02Cl2, in the framework of an extended one-band 
Hubbard model. By combining analytical and numeri- 
cal techniques, we found a consistent description of both 
charge and spin degrees of freedom responses as ob- 
served by optical conductivity, ARPES, magnetic Ra- 
man and inelastic neutron scattering experiments. We 
find that the almost dispersionless band measured by 
ARPES around (0, tt) ( relative to (7r/2,7r/2) ) on the 
one-hole dispersive mode may be ascribed to long-range 
hopping processes. We calculate also the non-resonant 
Big and Ai^ Raman profiles. The Big line is mainly 



composed of two nearby structures. One of them origi- 
nates on 2-magnon excitations and peaks at LU2mi while 
the other, centered around LUim, is due to 4-magnon scat- 
tering. For Sr2Cu02Cl2, we obtain uj2m ~ OSAeV and 
'-^im ~ O.GieV. The 4-magnon Raman signal, induced 
by multispin interaction terms, is a characteristic of the 
Hubbard model and has a negligible intensity for the min- 
imal two-dimensional AFH Hamiltonian. The Aig mode, 
by contrast to the AFH model, shows a finite Raman 
signal at frequencies around uj4m- To the best of our 
knowledge, this is the first time that an accurate calcu- 
lation of all those properties is obtained using the same 
model with the same parameter set. In section II, we 
describe the one-band Hubbard model used in this work 
along with the procedure that we follow to obtain the 
effective interaction parameters. Section HI and IV are 
devoted to the quasiparticle dispersion and to the spin 
excitations respectively, while in section V we summa- 
rize the results. 



II. THE EFFECTIVE MODEL: 

The effective one-band Hubbard model considered here 
includes the on-site Coulomb repulsion U and hoppings 
up to third-nearest neighbors. Microscopically, these 
hopping processes originate on the overlap between Wan- 
nier orbitals of a more complicated multiband model . 
Although, the strength of these interactions decreases 
with distance, recent work ||] suggests that hoppings fur- 
ther than first-nearest neighbors have to be included to 
obtain a quantitative description of experimental data 
for the cuprates. Of course, these hopping processes 
are material dependent. Here, we will focus our study 
on Sr2Cu02Cl2 lamellar cuprate, and then provide es- 
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timates for this material only. It is expected that al- 
though they introduce frustration and tend to decrease 
the strength of spin-spin correlations, their small values 
will not destroy the antiferromagnetic insulating ground 
state. However, they could play an important role on the 
charge dynamics. 

In standard notation, the dispersion for the kinetic en- 
ergy part of the single band Hubbard effective Hamilto- 
nian Cq is written as eq = V' + Cq +eq ^ +eq ' • Here, is a 

(r) 

constant and Cq , with r = 1,2,3, are the tight-binding 
dispersions for 1"*, 2"'' and 3'"^-nearest neighbors with 
hoppings — ti, t2 and respectively. For realistic values 
of multiband parameters, the effective hoppings t2 and 
t3 have the same order of magnitude of the corrections 
due to the states dropped by the reduction to a single 
band model Therefore, ti is the most appropriate 
energy scale. The value of ti was fixed at OAbeV 
while the other parameters were obtained by comparison 
with ARPES data. Our strategy is to solve first this dif- 
ficult many-body problem in a mean-field approximation 
and then, by using ARPES data, determine the value 
of t2, ts and the on-site Coulomb interaction U. Since 
Sr2Cu02Cl2 is an antiferromagnetic insulator, we use a 
spin-density wave (SDW) ansatz in the mean-field cal- 
culation. Notwithstanding its apparent simplicity, this 
treatment of the insulating half-filled t — U Hubbard 
model provides a successful description of the electronic 
degrees of freedom up to intermediate values of U . 
This analytical treatment of the Hubbard model has pro- 
vided also important inside in our current understanding 
of the resonant Raman scattering in antiferromagnetic 
insulators fl^ . 

The hole quasiparticle dispersion in the SDW approx- 
imation is given by. 



e(q) 



tq+TT 



(1) 



where n fixes the value of the Hubbard gap. Using 
ARPES data, we find n - 0.75ey, ^2 = 0.35ti, ^3 = 
O.OSii and -0 = Q.QQeV. The reduction of the three-band 
model onto the single band Hubbard model for realis- 
tic values of the parameters, indicates that the effective 
hopping ti is bounded between 0.3eV and 0.5eV while 
U /ti ~ 7—9. Furthermore, the derivation of the one band 
Hubbard Hamiltonian given by Simon and Aligia ( see 
Ref. 1^ ) for the parameters obtained from LDA calcula- 
tions for La2Cu04, gives ti ^ QAbeV, U/ti ~ 7.6, t2/ti ^ 
0.15, t^/ti ~ —0.12. Remarkably, the value of U obtained 
from the reduction agrees well with the one found from 
ARPES data. Note that differences in magnitude and/or 
sign between our estimates and the calculated t2 and 
are expected because they depend strongly on the sur- 
roundings of the Cu02 plane. 

Optical absorption measurements |^] on insulating 
Sr2Cu02Cl2 provide an additional check on k. These 



experiments show a charge-transfer absorption edge be- 
ginning at ~ 1.65eV and a strong band at w ~ 1.5eV. 
The latter was identified as an excitonic excitation. Re- 



cently, it was shown |13| that the observed absorption 
Eu peak lying at (0.1-0.2)eV below the absorption edge 
can be explained within an effective generalized one-band 
Hubbard model obtained from the simplest three-band 
model supplemented with the nearest-neighbor Coulomb 
interaction Upd- Aside from the on-site Coulomb in- 
teraction, this generalized Hubbard model includes the 
nearest- neighbor charge-charge interaction V . For sim- 
plicity, we have not taken into account either Upd nor V ^ 
and therefore the effective model consider in this work 
can not describe excitonic-like excitations. 

At the mean-field level, the optical conductivity Uxxi'^) 
does not depend on t2 and ts, and it is given, at T = 0, 
by 



. V 27r ^ ^ 9-27 ^ 
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where i?k 



,(Hl2 



The onset of the optical con- 



ductivity ffxxi^) found in the SDW approximation is at 
A = 2k ~ (1.50 ± 0.15)eT^, in agreement with the ex- 
perimental value for the charge-transfer absorption edge, 
A ~ 1.65eV. In the mean-field approximation, n is re- 
lated to the renormalized Coulomb parameter U through 
the mean- field gap equation ||ll|. For the ti — U Hub- 
bard model, one obtains U /ti = 1.80, 2.34, 5.80 for 
A/^i — 0.57, 1.05, 4.80. These values correspond to 
the hare Coulomb parameter U/ti — 2, 4, and 8 re- 
spectively. In all cases, U > U ||l^. In order to obtain 
the bare Coulomb interaction parameter [/, we calculated 
the gap by performing QMC simulations for different val- 
ues of J7, /3 = 5 — 12 and particle density < n >. By 
changing the doping from holes to electrons, the chem- 
ical potential fi crosses the gap at < n 1 where a 
plateau shows in the < n > vs /i curve, see Fig. la. On 
this plateau, the electronic compressibility K. vanishes, 
indicating an insulating state at that density. The width 
of the region with JC = measures the value of the charge 
gap and in turn allows us to provide an estimate of the 
bare Coulomb repulsion U. We found that the insulating 
Sr2Cu02Cl2 material can be described as an interme- 
diate coupling one-band Hubbard system with U/ti ~ 8 
and the other parameters as described above. Although, 
2"'' and 3'"'' nearest-neighbors hoppings introduce some 
degree of frustration on the magnetic background, the 
insulating ground state is still antiferromagnetic, as is 
found by performing QMC simulations on square clus- 
ters of (4x4) and (6x6) sites. In Fig. lb, we plot the mag- 
netic structure factor 5(q). The antiferromagnetic peak 
at (tt, tt) is clearly evident and its strength increases as 
the system size is increased, a signal of dominant antifer- 
romagnetic spin correlations. 
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FIG. 1. (a) < n > vs ^ at U/ti = 8. The flat region 
is a measure of tlie Mott-Hubbard gap. (b) Spin structure 
factor of tlie Hubbard model with hoppings up to third near- 
est-neighbors. QMC results for (full square) (4x4) and (open 
circle) (6x6) clusters. 

From QMC simulations, we also obtain the local mo- 
ment of the effective sites. Since an effective site rep- 
resents a Cu02 cell, an estimation of the local moment 
per Cu fj,s can be obtained by taking into account the 
Cu occupation on the cell. At stoichiometry, for a Cu 
occupation of ~ 80%, we obtain /ig 0.37/^s/Cu which 
is consistent with the experimental value jl^]. By con- 
trast, note that the measured local moment for insulating 
La2Cu04 is roughly twice the value of Sr2Cu02Cl2. 



III. QUASIPARTICLE DISPERSION: 

Soon after ARPES' results for the insulating cuprate 
Sr2Cu02Cl2, several theoretical works 16-1^ have been 
devoted to the description of the data by — J like hamil- 
tonians. Unfortunately, ARPES data show that the 2D 
ti~ J model accurately describes only e(q) along the di- 
rection from the zone center T = (0,0) to M = (7r,7r). 



Important differences were found moving along the non- 
interacting Fermi surface X = (tt, 0) — > (0, tt) and near 
the X point. 



Figure 2 shows a comparison of the single-hole disper- 
sion obtained from: ARPES data, ti — t^ — J model, and 
the Hubbard model with hoppings up to third nearest- 
neighbors. The theoretical e(q) can be obtained from an 
approximate treatment of the single hole problem. For 
the ti ^ t2 — J model, it can be determined by using 
the self-consistent Born approximation [ p^ . The quality 
of this approximation was contrasted successfully against 
ED calculations For the ti-t2-h-U model, the 

simplest procedure is to use the mean-field SDW analy- 
sis. Since, long-range hoppings are small in magnitude, 
we expect as for the ti~U model, the effect of quantum 
fluctuations can be absorbed into renormalized hopping 
values while the form of the dispersion relation remains 
the same as at the mean-field level. 



> 
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FIG. 2. Comparison between the quasiparticle dispersion 
of: the IBHM (solid line) with hopping up to third near- 
est-neighbors treated in the SDW approximation, ti — t2 ~ J 
(dashed line) and ti ~ J models (dot-dashed line) in the 
Born approximation for J = Q.125meV, J/ti — 0.3 and 
t2/ti = -0.35 |l|], and the ARPES data (full circles)obtained 
for insulating Sr2Cu02Cl2 §. 

The ARPES dispersion e(q) is described rather well by 
these theoretical models around the M point, possibly 
because they describe properly the magnetic structure 
of the quasiparticle cloud for this particular value of q. 
Outside the antiferromagnetic Brillouin zone, results for 
the ti —t2 — J model differ significantly from the experi- 
mental data, even for the case of finite and positive ^2/^1- 
This hopping process pushes e{X) down and at the same 
time decreases the bandwidth W |2^. Although sec- 
ond and further nearest-neighbor hoppings have a small 
strength, and at a first sight they seem to be irrelevant, 
they have important effects on the quasiparticle disper- 
sion, in particular around the X q-point and on the band- 
width's value. While a finite ^2 reduces the bandwidth, 
the main effect of is to increase W and, at the same 
time, reduce the dispersion around (0, tt). Meanwhile, 
without these interactions, it was found from QMC cal- 
culations that the t — U Hubbard model at U/ti = 8 
gives a bandwidth smaller than the experimental value 
(W = 280 meV) by a factor of 2. Including n.n.n. hop- 
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pings, we obtain an overall good description of the ex- 
perimental e(q) based on the functional form provided 
by the SDW mean-field solution. Although, we have not 
performed the highly demanding QMC computation of 
the single-hole dispersion for the model proposed in this 
work, our confidence on the SDW approximation comes 
from its success in describing the ti — U dispersion re- 
lation and the comparison performed in Sec. IV against 
other experiments using the very same parameter set. 
Further support is found from the recent calculation of 
the single hole dispersion done in Ref. ||2^ for the strong 
coupling limit of a generalized Hubbard model. Note 
nevertheless, that in this calculation a J bigger {^-^ 17% 
) than the experimental exchange constant is required to 
obtain the experimental band width. Of course, within 
the SDW scheme we can not study other interesting prop- 
erties of the quasiparticle such as its residue |^ . 

Along the non-interacting Fermi surface, results ob- 
tained from Hubbard like models are in better agreement 
with ARPES measurements than the one hole ti — t2 — J 
dispersion. Let us emphasize that only a few experimen- 
tal points, taken from panel (a), were used to determine 
the hopping parameters. As a by-product, the theoret- 
ical dispersion agrees also rather well with ARPES re- 
sults of panel (b) . The small asymmetry observed along 
the (0, tt) (tt, 0) line could be ascribed to sample 
anisotropics. 

IV. SPIN EXCITATIONS: 

On the experimental side, the spin degrees of freedom 
are tested by Raman and neutron scattering experiments. 
They reveal, in fact, that the insulating ground state 
of Sr2Cu02Cl2 is antiferromagnetic. The experimental 
value of the spin wave velocity is c ~ 0.83 (el^ — A) 
. At low temperatures -g, spin waves excitations con- 
tribute to the internal energy per site e(/3). Following 
Tang and Hirsch (25|, we first calculate e(/9) using the 
QMC method and then by fitting the spin- wave contribu- 
tion to the internal energy, we estimate c. While for the 
t—U Hubbard model, the spin wave velocity for U/ti = 8 
(c ~ 1.10(ey — A) is bigger than the experimental value, 
for the parameter set proposed for Sr2Cu02Cl2 we find 
c ~ 0.85{eV — A) in fairly good agreement with the avail- 
able data. 

The scattering of light from insulating antiferromag- 
nets at a low energy scale compared with the charge- 
transfer gap A, provides additional information about 
the spin dynamics. The shape of the Big Raman profile 
R{lu) has interesting features, namely, a characteristic 
peak ascribed to 2-magnon excitations, a broad linewidth 
and a very asymmetric profile with a " shoulder-like" fea- 
ture at higher frequencies, but close to the 2-magnon 
peak. At a first sight, the two latter features seem to 
be mainly due to different physical phenomena, namely 



spin-phonon interaction |2^] and quantum spin fluctu- 
ations respectively |2^]. Evidence for other contribut- 
ing mechanisms to the width of the 2-magnon line, aside 
from the quantum spin fluctuations, comes from the fact 
that the half-width of the Big Raman response has al- 
most the same value ^ 1200 cm^^ for all members of the 
M2CUO4 series although the exchange constant changes 
by ~ 20%, i.e. the width of the 2 — magnon line does not 
scale with J. Furthermore, it was argued recently that 
the spin-phonon interaction can be responsible for the 
broad linewidth observed on this geometry ||2^. In fact, 
ED and QMC calculations of the Raman cross section 
on the 2D-AFH model supplemented with spin-phonon 
interactions describe the broad linewidth observed in the 
insulating compounds of high-Tj. superconductors. De- 
spite the theoretical success in describing the position 
and linewidth of the 2-magnon line, current results sug- 
gest that the description of the " shoulder- like" feature, 
whose position was assigned experimentally to ~ 4J, 
require to go beyond the minimal AFH model. In fact, 
a detailed study of the effect of four-magnon scattering 
in the 2D- AFH model shows that the intensity of the 
Raman signal results too small to fully account for the 
experimental data pOt . 

Additional terms ( multispin interactions) appear quite 
naturally from the one-band Hubbard model scheme. In 
fact, by performing a canonical transformation up to 4*'' 
order on the Hubbard model one obtains an effec- 
tive spin Hamiltonian which besides the antiferromag- 
netic exchange interactions up to third-nearest neigh- 
bors, includes a 4-spin cyclic exchange term with strength 
^ 80tf/U^. At U/ti = 8, the exact and effective ground 
state energies differ by less than 1%( see S.Bacci et al. in 
Ref. 1^). In Fig. (3), we plot the non- resonant Big Ra- 
man spectrum obtained from ED calculations on a ^20 
X \/20 cluster. In this calculation, we use the traditional 
Hamiltonian for describing the interaction of light with 
spin degrees of freedom, i.e. the Loudon-Fleury Hamilto- 
nian, which in standard notation is written as 

Obi = ^ <5'i.(S'i+e^ - S'i+ey). (3) 
i 

and the now standard continued fraction approach p^ ] 
to obtain the Raman line. Although, we did not perform 
flnite-size scaling, finite-size effects are small because of 
the local nature of the Raman operator. The calculation 
of the resonant scattering contribution to the Raman sig- 
nal is out to the scope of this work. As for the resonant 
case iQ, the non-resonant Big profile is composed of two 
structures, namely a 2-magnon peak at uJ2m ~ O.^AeV 
and a side band centered around ~ 0.64ey, in rea- 
sonable agreement with the experimental values . The 
Raman signal around is mainly due to 4-spin cyclic 
exchange interaction terms. 
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V. SUMMARY 
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FIG. 3. The Big (a) and Aig (b) non-resonant Raman spec- 
tra of the efTective spin Hamiltonian. The value of the param- 
eters are as in Fig.(l). (a) Dashed lines are for the experi- 
mental results of Ref. and the solid line is for the 20 sites 
cluster, (b) Aig line for clusters of 16, 18 and 20 sites. 

The first moment of this line is Mi ~ OAeV. Within 
the context of the AFH model, we obtain from Mi 
{— 3.6 Je ||]), the exchange constant Jg ~ lllmeV which 
is roughly consistent with Jg ^ l25meV as inferred from 
neutron scattering [||. For the Hubbard model, the 2- 
magnon excitation energy depends not only on the bare 
exchange constant J ^ Atf/U — 2At\/U^ but also on the 
degree of frustation introduced by 2"'' and 3'"'' neighbors 
exchange processes. These terms produce a shift of the 
peak towards zero frequency |^ and as a consequence 
a strong renormalization of the microscopic J could take 
place. Our results, based on the Hubbard model, sug- 
gest that J is almost twice the effective Je- For the Aig 
symmetry, the Raman operator given by 

O^, -^S'i.(5i+e.+S'i+eJ. (4) 
i 

does not commute with the effective spin Hamiltonian 
and produces a finite signal in this otherwise forbidden 
channel. The Aig line shape is very asymmetric with al- 
most all the spectral weight around Wim- At higher fre- 
quencies, multimagnon scattering gains intensity, making 
this line broader than R{oj) for the Big symmetry. 



In summary, our observations and conclusions sup- 
port previous analytical work base on a systematic low- 
energy reduction of complicated multiband onto a single- 
band Hubbard model. We find that a single-band Hub- 
bard model supplemented with hoppings up to 3'"'^ near- 
est neighbors describes several experimental features ob- 
served on insulating Sr2Cu02Cl2. Let us emphasize that, 
while the parameters of this single-band model were de- 
termined from a few experimental ARPES data points 
( not a fit ), i.e. charge degrees of freedom, we were also 
able to describe as well spin excitations. Our results for 
the quasiparticle dispersion resemble the ARPES disper- 
sion and suggest that the almost dispersionless part ( 
relative to (7r/2,7r/2) ) around (0, tt) could be ascribed 
to long-range hopping processes. Although, t2 and ^3 
introduce frustration on the magnetic background, the 
system is still an antiferromagnet. Our results for the 
description of spin excitations have implications for the 
interpretation of the mid-infrared optical absorption in 
undoped lamellar copper oxides Multispin terms, 
introduce multimagnon processes that, could contribute 
significantly to the weight of the sidebands |^ . 
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